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Magic state distillation is an important primitive in fault-tolerant quantum computation. The magic states 
are pure non- stabilizer states which can be distilled from certain mixed non- stabilizer states via Clifford group 
operations alone. Because of the Gottesman-Knill theorem, mixtures of Pauli eigenstates are not expected to be 
magic state distillable, but it has been an open question whether all mixed states outside this set may be distilled. 
In this Letter we show that, when resources are finitely limited, non-distillable states exist outside the stabilizer 
octahedron. In analogy with the bound entangled states, which arise in entanglement theory, we call such states 
bound states for magic state distillation. 
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The significant noise and decoherence in quantum systems 
means that harnessing these systems for computational tasks 
must be performed fault tolerantly |1, 2|. In a wide variety 
of setups only a limited set of gates, known as the Clifford 
group, are implemented in a manifestly fault tolerant man- 
ner. Examples include some anyonic topological quantum 
computers ||3}Q, post-selected quantum computers |[6j[3 and 
measurement based topological quantum computers fSl. This 
motivates the problem of when such devices, with practically 
error free Clifford gates, may be promoted to a full quantum 
computer. The celebrated Gottesman-Knill theorem shows 
that a Clifford circuit acting on stabilizer states — simultane- 
ous eigenstates of several Pauli operators — can be efficiently 
simulated by a classical computer ||9l. However, given a re- 
source of pure non-stabilizer states, we can implement gates 
outside the Clifford group. For example, a qubit in an eigen- 
state of the Hadamard enables one to implement a 7r/8 phase 
gate that when supplementing the Clifford group gives a dense 
covering of all unitary operations [ lO L and so enables univer- 
sal quantum computation. 

Preparation of non- stabilizer states would usually require a 
non-Clifford operation, so in this context, one would require 
that even noisy copies of these states enable high fidelity quan- 
tum computation. Bravyi and Kitaev fTOl showed that this can 
be achieved. Coining the term magic state distillation, they 
showed that most mixed non- stabilizer states can be distilled 
via Clifford group circuits to fewer copies of a lower entropy 
state, reaching in the limit of infinite iterations a pure non- 
stabilizer magic state. However, the protocols they presented 
do not succeed for all mixed non- stabilizer states. Bravyi and 
Kitaev were not satisfied by the ambiguous status of these 
states and concluded that ''The most exciting open problem is 
to understand the computational power of the model in [this] 
region of parameters.''. Either all non- stabilizer states are ef- 
ficiently distillable by an undiscovered protocol, or there exist 
non- stabilizer states that are impossible to distill. Such undis- 
tillable states we call hound states for magic state distillation, 
in analogy with bound states in entanglement distillation ifTTl . 
Here we make progress by showing that bound states exist 
for a very broad class of protocols. By showing that a single 
round of a finite sized protocol will not improve these states, 
it follows that repeating such a protocol, even with an infinite 



number of iterations, will also have no benefit. Hence, we 
explain why all known protocols fail to distill some states. 

The single-qubit stabilizer states, for which the Gottesman- 
Knill theorem applies, are the six pure stabilizer states (the 
eigenstates of ±X, ±F and ±Z) and any incoherent mixture 
of these. In the Bloch sphere, this convex set with 6 vertices 
forms the stabilizer octahedron partially shown in figure [T^. 
Single-qubit states have density matrices: 

a) = (1 + (2/ - l){axX + ayY + azZ)) /2, (1) 

where a = (ax, <^y , «z) is a unit vector, and / is the fidelity 
w.r.tthepurestate |V^a)(^a| = {t ^ axX ^- ayY ^ azZ) /2. 
Stabilizer states satisfy: 

|2/-l|(|ax| + |ay| + |az|)<l (2) 

where the equality holds for states on the surface of the oc- 
tahedron, and we denote the fidelity of such surface states as 
f^, which is unique assuming / > 1/2. 

Prior protocols for magic state distillation fJ] [TOl [121 [13 
increase fidelity towards eigenstates of Clifford gates, such as 
the Hadamard H and the T gate |[T8ll . These eigenstates have 
SiH = (1,0,1)/V2 andaT = (1,1,1)/V3, with / = 1 for 
ideal magic states. Given the ability to prepare a mixed non- 
stabilizer state, p, we can perform an operation called polar- 
ization, or twirling, that brings p onto a symmetry axis of the 
octahedron. For example, by randomly applying 1, T or T^, 
we map p p{f^ a^). 

Bravyi and Kitaev proposed the following protocol LI 01 for 
|T) state distillation: (1) Prepare 5 copies of p(/, a^); (2) 
Measure the 4 stabilizers of the five-qubit error correcting 
code; (3) If all measurements give +1, the protocol succeeds 
and the encoded state is decoded into a single qubit state, and 
otherwise restart. Upon a successful implementation of this 
protocol the output qubit has a fidelity F{f) plotted in fig- 
ure[2]3. Provided the initial fidelity is greater than some thresh- 
old, a successful implementation yields a higher fidelity. This 
protocol has a non-tight threshold, and exhibits a gap between 
the threshold and the set of stabilizer states. Because the initial 
state was twirled onto the T axis, the threshold forms a plane 
in the Bloch sphere (see figure [T]). In contrast, Reichardt has 
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FIG. 1: One octant of the Bloch sphere with various regions and 
directions shown, (a) The blue region shows the stabilizer states in 
one octant. Each octant is identical, with all stabilizer states forming 
an octahedron, (b) The yellow plane is the distillation threshold for 
the 5 qubit code, with the direction of a.T shown. The yellow plane 
is parallel to the underlying blue face of the stabilizer octahedron, 
but is displaced by a small gap. (c) The three green planes are the 
thresholds for the Steane code, with each plane differing by local 
Clifford gates. These planes meet the stabilizer octahedron at its 
edges. The three vectors are axes of H-like gates, e.g. slh', (d) The 
combined region of states distilled by either the 5 qubit code or the 
Steane code. This region only touches the stabilizer octahedron at its 
edges, and no other known protocol is tight in any other direction. 



proposed a protocol that does have a tight threshold for distil- 
lation of S-h) states in a i^-like direction 1 12 1. His pro- 
tocol is similar to above, but uses 7 qubits each attempt and 
measures the 6 stabilizers of the STEANE code |1|. In fig- 
ure [2^ we show the performance of this protocol, where there 
is no threshold gap. When the initial mixture is not of the form 
S-h), we twirl the initial mixture onto the H axis. Hence, 
the threshold forms a plane for each i^-like direction (see fig- 
ure [T]). Although the protocol is tight in directions crossing 
an octahedron edge, the protocol fails to distill some mixed 
states just above the octahedron faces, and so is not tight in all 
directions. Even the combined region of states distilled by all 
known protocols still leaves a set of states above the octahe- 
dron faces, whose distillability properties are unknown. 

Here we show that for all size n protocols there is a re- 
gion of bound states above the octahedron faces. More for- 
mally, we considering all states ap) where ap has all 
positive (non-zero) components. Having all components as 
non-zero excludes states above octahedron edges. Consider- 
ing only states in the positive octant is completely general 
as Clifford gates enable movement between octants. Many 
copies of bound states cannot be used to improve on a single 
copy, and below we formalize the idea of not improved and 
state our main result. 

Definition 1 We say p' is not an improvement on p(/, ap), 
when p' is a convex mixture of Cip{f^ap)Cl and stabilizer 
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FIG. 2: The performance of magic state distillation of: (a) the 
STEANE code for distilling states in a if-like direction; (b) the five 
qubit code distilling states in a T-like direction. Notice that both 
functions are continuous, and that in (b) an input state on the oc- 
tahedron surface, / = /a^, will output a state below the surface. 
Consequently, there is also a region above the surface where the out- 
put is a stabilizer state, and hence, not an improvement on the initial 
state. 



states, where Ci are Clifford group gates. 

Theorem 1 Consider a device capable of ideal Clifford gates, 
preparation of stabilizer states, classical feedforward and 
Pauli measurements. For any protocol on this device that 
takes p{f,8ip)^^ and outputs a single qubit, p', there exists 
an e > such that p' is not an improvement on p{f^ ap) for 
f < ft + e. 



Theorem [T] covers a wide class of protocols, which attain 
a fidelity that is upper-bounded by a narrower class of proto- 
cols |[T4iL such that theorem [T] follows from: 

Theorem 2 Consider all protocols that follow these steps: ( i) 
prepare p{f^ ap)*^^; (ii) measure the n — 1 generators of an 
n qubit stabilizer code Sn-i ^ith one logical qubit; (Hi) post- 
select on all '*+7" measurement outcomes; (iv) decode the 
stabilizer code and output the logical qubit as the single qubit 
state p' . For all such protocols there exists an e > such that 
p' is not an improvement on p{f^ ap)for f < + e. 



Prior protocols, such as those based on the STEANE code 
and 5 qubit code, are covered explicitly by theorem [2] Here 
we use the structure of stabilizer codes to prove theorein[2j 
with theorem [T] following directly from the results of (T4J, 
where such distillation protocols are shown to have equal ef- 
ficacy with more general Clifford protocols. It is crucial to 
consider the implication of these theorems when an n-qubit 
protocol is iterated m times. When a single round provides no 
improvement on the initial resource, the input into the second 
round will only differ by Clifford group operations, and hence 
our theorem applies to the second, and all subsequent, rounds. 
Hence, repeated iteration cannot be used to circumvent our 
theorem. Before proving these theorems, we derive a pair of 
powerful lemmas that identify bound states. 

Lemma 1 Consider n copies of an octahedron surface state 
p{f^^,ap) projected onto the codespace of Sn-i cind then 
decoded. If the output qubit is in the octahedron interior, then 
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there exists an e > such that for f < f^^ + e the same 
projection on p{f, s.p)'^^ also projects onto a mixed stabilizer 
state. 

This lemma follows directly from the dependence of the out- 
put on /, which for finite n is always continuous. We can 
observe this lemma at work in figure |2]). Our next lemma 
identifies when octahedron surface states are projected into 
the octahedron interior. Before stating this we must establish 
some notation. An initial state, p{f^^ , ap)*^"^, is an ensemble 
of pure stabilizer states: 

/>(/fp,ap)^-= ^ ^g|^g)(^g|, (3) 

gG{x,y,z}- 

where |^g) is stabilized, ^|^g) = |^g), by the group Qg 
generated by g =(gi, g2,--- Qn)- The operator Qi is X^, Yi 
or Zi, with i labeling the qubit on which it acts. Each con- 
tribution has a weighting = Yli{cig-/{ax + cly + clz))- 
Measuring the generators of Sn-i and post-selecting on 
outcomes, projects onto the codespace of 5^-1 with projector 
P = T^seSr^-i producing: 

with projected terms, |^g), of new weighting q'^. Each |^g) 
has its stabilizer generated by (Gg, si, S2, ....5^-1), where Gg 
is an independent generator that: (a) was present in the ini- 
tial group Gg G ^g; and (b) commutes with the measure- 
ment stabilizers GgSn-i = Sn-iG^. In other words, it 
must be equivalent to one of six logical Pauli operators of 
the codespace. We denote the set of logical operators as £, 
and its elements iX^, ±Yl and iZ^, and so Gg G C.Sn-i- 
This defines a decoding via the Clifford map, Xl ^ Xi and 
Zl ^ Zi. Since there are only six distinct logical states, we 
can combine many terms in equation |4j 

where \^ l) has stabilizer generators (L, 5i, S2, ....5^-1). The 
new weighting \^ qL = with the sum taken over all g 

that generate Qg containing an element Gg G L.Sn-i- We 
can now state the next lemma: 

Lemma 2 Given n copies of p{f^^,ap) projected into the 
codespace of Sn-i cind decoded, the output qubit is in the 
octahedron interior if there exist any two pure states in the 
initial ensemble, |^g) and |^g/) (defined in equation^, such 
that both: 

( i) the projected pure states are orthogonal, so that L ^ Qg 
and —sLe Q^/ where L ^ C and s G Sn-i; 

and 

(ii) upon projection |^g) and |^g') do not vanish, so q'^ ^ 
and q' f ^ 0. 



We prove this lemma by contradiction. From equation [2j and 

(2/ — l)aL = {qL — Q-l), surface states satisfy: 

\qxL - Q-Xl I + \QYl - Q-Yl I + \QZl - Q-z^ I = 1, (6) 

and we assume to the contrary that the projected state has this 
form. Since q±L are non-negative reals, we have \qL — Q-L \ = 
Ql -\-Q-l — 2Min(g'L, ^-l), where Min(g'L, ^-l) is the mini- 
mum of qL and q-L- Along with the normalization condition, 
^i^qL = 1, this entails: 

Mm{qxL , q-XL ) + Min((7y^ , q_Y^ ) + Miniqz^ ,q_ZL) = ^- 

Since all terms are positive, no cancellations can occur and 
so every term must vanish, hence Min(gL, g-^) = 0,VL. 
However, conditions (/) and (//) of the lemma entail that there 
exists a non-vanishing Min(g'L, g-^), as > ^g 7^ and 
q-L > ^g/ 7^ 0. Having arrived at this contradiction, we 
conclude the falsity of the assumption that the projected state 
remains on the octahedron surface, and so must be in the octa- 
hedron interior. This proves lemma [2j and we now show that 
lemma[2]applies to all stabilizer reductions that do not trivially 
take/,(/,ap)®"^ap(/,ap)C|. 

Our proof continues by finding canonical generators for the 
code Sn-i' A related method has been used to prove that 
all stabilizer states are local Clifford equivalent to a graph 
state 1 15], and we review this first. All stabilizer states have 
a stabilizer Sn with n generators. Each generator is a tensor 
product of n single-qubit Pauli operators. This can be visual- 
ized as an n by n matrix with elements that are Pauli opera- 
tors, each row a generator and each column a qubit. Different, 
yet equivalent, generators are produced by row multiplication, 
via which we can produce a canonical form. In this form col- 
umn i has a non-trivial Pauli operator Ai that appears on the 
diagonal, and all other operators in that column are either the 
identity or another operator Bi. Note that Ai and Bi compose 
a third non-trivial Pauli AiBi = i{—l)^^Gi with 7^ = 0, 1. 
Hence, all stabilizer states differ from some graph state by 
only local Cliffords that map (A^, Bi) (X^, Zi). 

A code, Sn-i, has one less generator than the number of 
qubits, and so more columns than rows. We can apply the 
diagonalisation procedure on an n — 1 by n — 1 submatrix, to 
bring this submatrix into canonical form. Hence, we can find 
generators of 5n-i such that: 

s,={-ir^A,{l[ Bt')T,,n, (7) 

where the variables ^kj = 0, 1 denote whether B^ or tk is 
present, and aj =0,1 defines the phase. With the n*^ column 
out of canonical form, this leaves the n*^ qubit operator Tj^n 
unspecified. However, if all these generators have Tj^n = ^n. 
then the protocol is trivial and projects n — 1 qubits into a 
known stabilizer state and the last qubit untouched, and so no 
improvement is made for any /. Hence, herein we assume the 
non-trivial case; in particular we assume stabilizer Tn-i^n 7^ 
In. Since, we can always relabel qubits this is completely 
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general. Furthermore, we can define T^-i,n = ^n- Now we 
can define a logical operator in the codespace of 5^-1 : 

= I n Bn-lBn, (8) 

\l<j<n-2 J 

where the variables Q = 0, 1 are uniquely fixed by commu- 
tation relations ZlSj = sjZl. Note that Zl has some inbuilt 
freedom as Bn is not fixed other than that Bn 7^ A^, In, 
which is equivalent to free choice of 7^ in the expression 
AnBn = z( — Now wc enquire whether the final state 
contains two terms stabilized by and —sZl respectively, 
hence satisfying the conditions for lemma [2] If we consider 
the product of Zl and Sn-i, and choose 7^ = o^n-i + 7n-i 
mod 2, we have: 

- Sn-lZL = I n ^f''-^^^' I Cn-lCn. (9) 
yi<fe<n-2 J 

Our choice of 7^ ensures a minus sign on the left hand side, 
which aids in finding |^g) and |^g') that satisfy our lemma 
by being stabilized by Gg = Zl and Gg/ = —Su-iZl re- 
spectively. This criterion is fulfilled when: 

g = {Bi,B2, , Bn-2, Bn-l, Bn), 

= (^1,^2, , ^n-2, Gn-l, Gn). 

These states only vanish under projection, g'g, q^^, = 0, if they 
are stabilized by the negative of some element of the code 
Sn-i- To prove they don't vanish, we first observe that every 
element of Qg and ^g/ has either tj or Bj acting on qubit j, 
for all j = 1, 2, ...n — 2. The only elements of 5n-i for which 
this is true are t and Sn-i, but Sn-i has A^-xA^ acting on 
the last two qubits and neither or Q^i contain any such 
element. 

Using a canonical form of the generators of 5^-1, we have 
shown that non-trivial codes always satisfy the conditions of 
lemma [2] That is, all non-trivial codespace projections take 
many surface states into the octahedron interior. From the 
continuity expressed by lemma [T] this entails the existance of 
a finite region of non- stabilizer states that are also projected 
into the octahedron. Hence, all n-copy protocols do no im- 
prove on a single copy for some region of bound states above 
the octahedron faces, completing the proof. This does not con- 
tradict known tight thresholds in edge directions, as these di- 
rections have a with one zero component. 

Although our proof holds for protocols using fixed and fi- 
nite n copies of ap), we could conceive of a protocol that 
varies n. If this varying-n protocol has an n-dependent thresh- 
old, /J^(n), and /J^(n) f^^ as n 00, then its thresh- 
old would be arbitrarily suppressible. Repeated iterations of a 
protocol, or equivalently employing concatenation of a single- 
qubit code, will not change the threshold. However, one could 
consider a broader class of protocols consisting of iterates that 
act on p qubits and output q qubits (for p > q > 1) followed 



by a final round outputting a single qubit. Such protocols map 
n qubits to 1 qubit, with n growing each iterate, but with only 
p qubits involved in each iterate. This implies that multi-qubit 
output iterates may suppress the threshold effectively, and are 
worth further study. Currently, no such protocol is known. 
As such, in the asymptotic regime, bound magic states may 
not exist. However, numerical evidence so far indicates that 
smaller codes tend to produce better thresholds than larger 
codes. Nevertheless, the theorem does not rule out infinite 
cases from attaining a tight threshold. In the regime of finite 
resources, bound states do exist, and it is interesting ask what 
computational power Clifford circuits acting on such states 
possess. Can we find methods of efficiently classically simu- 
lating bound states; or can bound states be exploited in algo- 
rithms that offer a speedup over classical computation? 

Furthermore, our proof assumes a protocol acting on identi- 
cal copies, which invites study into whether our results extend 
to non-identical copies. In particular, following the analogy 
with entanglement distillation, we speculate that bound magic 
states may be distillable via "catalysis", where some non- 
consumed distillable resource activates the distillation ifTTIl . 
Finally we note that noisy Clifford gates can also enable quan- 
tum computation [16, 17 1, and we conjecture that a similar 
theorem will apply to a class of noisy Clifford gates analo- 
gous to states just above the octahedron faces. 
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